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We address an important issue as to whether bulk-sensitive data of Raman scattering, optical con- 
ductivity, magnetic penetration depth, directional point-contact tunneling spectra, and nonmagnetic 
pair-breaking effect in optimally electron-doped Ndi.ssCeo.isCuCU-;, support a nodeless s-wave or 
d-wave superconducting gap. We numerically calculate Raman intensities, directional point-contact 
tunneling spectra, and nonmagnetic pair-breaking effect in terms of both s-wave and d-wave gap 
symmetries. We find that all these bulk-sensitive data are in quantitative agreement with a nearly 
isotropic s-wave gap. The fact that T c is nearly independent of the residual resistivity rules out any 
d-wave gap symmetry. 
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The gap symmetry of high-temperature cuprate super- 
conductors has been a topic of intense debate for over 
twenty years. For hole-doped cuprates, bulk-sensitive 
experiments probing low-energy excitations in the su- 
perconducting state have consistently pointed towards 
the existence of line nodes in the gap function of hole- 
doped cuprates [H-Q, which is consistent with either d- 
wave gap (having four line nodes) or extended s-wave 
gap (having eight line nodes). Other bulk-sensitive ex- 
periments on hole-doped cuprates 0, Q can be quanti- 
tatively explained by extended s-wave gap 0j- In con- 
trast, the issue as to whether there exist line nodes in the 
gap function of electron-doped (n-type) cuprates remains 
controversial. Phase and surface-sensitive experiments 
Q provided evidence for pure d-wave order-parameter 
(OP) symmetry in optimally doped and overdoped n- 
type cuprates. Surface-sensitive angle-resolved photoe- 
mission spectroscopy (ARPES) @, [l(| showed a d-wave 
gap with a maximum gap size of about 2.5 meV. This 
gap size would imply a T c of about 14 K at the sur- 
face, which is significantly lower than the bulk T c of 
26 K (Ref. (if]]). Earlier magnetic penetration depth 
data [llj of optimally electron-doped cuprates with T c 
— 24 K were shown to be consistent with nodeless s- 
wave gap symmetry. Later on, a T 2 dependence of the 
penetration depth at low temperatures was observed in 
Pri. 85 Ceo.i5Cu04_ y (T c — 20 K), which appears to sup- 
port d-wave gap symmetry in the dirty limit 12]. But 
the same data can be quantitatively explained [7] in terms 
of a nodeless s-wave gap if one takes into account an ex- 
trinsic effect due to current-induced nucleation of vortex- 
antivortex pairs at defects. Extensive penetration depth 
data [l3j of Pr2- :c Ce :E Cu04-i, confirmed the nodeless gap 
symmetry at all the doping levels except for a deeply un- 
derdoped Pri.885Ceo.ii5Cu04_ y sample with T c = 12 K. 
Point-contact tunneling spectra (l4l- 18| also showed no 
zero-bias conductance peak (ZBCP) at all the doping lev- 
els except for a deeply underdoped Pri.87Ceo.i3Cu04_ y 
with T c = 12 K. Therefore, the penetration depth and 
point-contact tunneling spectra consistently suggest that 
the gap symmetry in deeply underdoped samples should 



be d-wave and change to a nodeless s-wave when the 
doping level is above a critical value. This scenario can 
naturally explain the d-wave gap symmetry inferred from 
surface-sensitive experiments if surfaces or interfaces are 
deeply underdoped. Experiments on hole-doped cuprates 
III 20] indeed show that surfaces and interfaces are sig- 
nificantly underdoped. 

Here we focus on optimally electron-doped 
Ndi.85Ceo.i5Cu04_j / (NCCO) to unambiguously address 
this important issue as to whether bulk-sensitive data 
of Raman scattering, optical conductivity, magnetic 
penetration depth, directional point-contact tunneling 
spectra, and nonmagnetic pair-breaking effect support a 
nodeless s-wave gap. We numerically calculate Raman 
intensities, directional point-contact tunneling spectra, 
and nonmagnetic pair-breaking effect in terms of both 
s-wave and d-wave gap symmetries. We find that all 
these bulk-sensitive data are in quantitative agreement 
with a nearly isotropic s-wave gap. The fact that T c is 
nearly independent of the residual resistivity rules out 
any d-wave gap symmetry. 

Bulk-sensitive Raman scattering has been proved to be 
a very powerful tool to study the anisotropy of the su- 
perconducting energy gap. Experiments carried out with 
different polarization orientations pick up the contribu- 
tions to the light scattering on different parts of the Fermi 
surface. The B\ g spectra provide information on the light 
scattering primarily in the neighborhood of the k x and k y 
axes while Bi g spectra probe mainly along the diagonals, 
where {k x , k y ) = k is the in-plane wave vector of elec- 
trons. The electronic Raman intensity is proportional to 
the imaginary part of Raman susceptibility X 77 (<f, w) in 
the limit of q approaching 0. At zero temperature, the 
imaginary part of x 77 (cj) is given by [2l| 
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where A(fc) is the momentum-dependent superconduct- 
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FIG. 1: Raman intensities at 8 K for Ndi.ssCeo.isCuCU-y (T c 
= 22±1 K) in the B\ g (Fig. la) and j?ga (Fig. lb) symmetries. 
The data are digitized from Ref. [24|. The solid lines are 
numerically calculated curves in terms of an anisotropic s- 
wave gap function: A = 3.59(1 — 0.11 cos AO) meV with V 
— 1.4 meV for the B\ g channel and V — 1.9 meV for the 
£?2 9 channel; the dashed lines are the numerically calculated 
curve in terms of a non-monotonic d-wave gap function: A = 
3.25(1.43 cos 20- 0.43 cos 60) meV with T = 1.6 meV for both 
channels. 



ing energy gap, E(k) = y A 2 (fc) + e 2 (k), T is the param- 
eter associated with the life-time broadening of the quasi- 
particles, e(k) is the band-dispersion relation, and 7(fc) 
is the Raman vertex which is proportional to COS k x Q> — 
cos k y a for B\ g symmetry and to sin k x a s\a.k v a for B^ g 
symmetry (where a is the lattice constant) [22J. By fit- 
ting ARPES data, the band-dispersion relation e(k) (in 
units of meV) for Ndi ssCeo lsCuC^-y was found to be 
231 



cos k y a) 



220 cos k x a cos k y a 



e(k) = — 460(cosfc x a - 

— 70(cos 2k x a + cos 2k y a) + 60(cos 2k x a cos k y a 

+ cos k x a cos 2k y a) — 27 (2) 



and B~2 g Raman intensities at 8 K 
with T c = 22±1 K. The data are 
digitized from Ref. [24]. The solid lines are numerically 
calculated curves using Eqs. 1 and 2, T = 1.4 meV for the 



Figure 1 shows B\ g 
for Ndi. 85 Ceo.i 5 Cu04- 



B\ g channel, F = 1.9 meV for the Bi g channel, and A = 
3.59(1 — 0.11 cos46>) meV, where 8 is measured from the 
Cu-0 bonding direction. It is remarkable that the calcu- 
lated curves match very well with the experimental data. 
Furthermore, the minimum value A m i„ of this gap func- 
tion is 3.2 meV, which is close to that (3 meV) deduced 
from the magnetic penetration depth [111 ]. The optical 
reflectance of a similar Ndi.ssCeo.isCuO^^ crystal with 
T c = 23 K also shows a minimum gap of about 3.1 meV 
at 10 K (Ref. [25[). For comparison, we also numerically 
calculate Raman intensities in terms of a non-monotonic 
ci-wave gap function: A = 3.25(1.43 008 26* — 0.43 cos 6(9) 
meV (dashed lines). This gap size, which is a factor of 
1.7 larger than that extracted from ARPES matches 
the peak position of the Bi g Raman spectrum. It is ap- 
parent that the Bi g Raman spectrum is inconsistent with 
any ci-wave gap function with nodes along the Cu-Cu di- 
rections. 

Further evidence for the nodeless s-wave gap symme- 
try comes from the directional point-contact tunneling 
spectra of optimally doped Ndi.ssCeo.isCuOzi-j, (T c = 
25 K). It was argued that single-particle tunneling ex- 
periments along the CuC>2 planes can probe the bulk 
electronic density of states since the mean free path is far 
larger than the thickness of the possibly degraded surface 
layer 26| . Figure 2 shows the in-plane point-contact tun- 
neling spectra of Ndi.ssCeo.isCuC^-j, measured along 
(100) direction (Fig. 2a) and (110) direction (Fig. 2b), 
respectively. The data are digitized from Ref. 17]. One 
can calculate point-contact tunneling spectra using the 
Blonder-Tinkham-Klapwijk (BTK) theory (27|. In this 
model, two parameters are introduced to describe the 
effective potential barrier (Z) and the superconducting 
energy gap A. As a supplement, the quasiparticle en- 
ergy E is replaced by E — iT, where T is the broadening 
parameter characterizing the finite lifetime of the quasi- 
particles. Based on the BTK theory, Shan et. al. calcu- 
late the tunneling conductance in terms of the isotropic 
s-wave gap function (l7| . The agreement between the 
calculated curve and data is excellent for each tunneling 
spectrum. However, the gap sizes that used to fit the 
tunneling spectra along the (100) and (110) directions 
are slightly different. This implies that the gap is not 
isotropic. 

For the extended anisotropic BTK model 28] , another 
parameter a is introduced to distinguish between differ- 
ent tunneling directions. In Fig. 2, we compare the tun- 
neling spectra with the calculated curves based on the 
extended anisotropic BTK model and an anisotropic s- 
wave gap function: A = 3.52(1 - 0.17 cos 40) meV. The 
solid line in Fig. 2a is the numerically calculated curve 
using Z = 2.8, T = 0.86 meV, and a = and the solid line 
in Fig. 2b is the calculated curve using Z = 3.0, T — 0.65 
meV, and a — 7r/4. It is apparent that the calculated 
curves are in excellent agreement with the data. 

In Figure 3, we compare the tunneling spectra with 
the calculated curves in terms of a nonmonotonic d-wave 
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FIG. 2: In-plane point-contact tunneling spectra of 
Ndi.85Ceo.isCu04-» ( T c = 25 K) measured along (100) di- 
rection (Fig. 2a) and (110) direction (Fig. 2b), respectively. 
The data are digitized from Ref. fl7| . The solid line in Fig. 2a 
is the numerically calculated curve using Z = 2.8, F = 0.86 
meV, and a — and the solid line in Fig. 2b is the calculated 
curve using Z = 3.0, T — 0.65 meV, and a = n/4. The gap 
function used in the calculations is A = 3.52(1 — O.17cos40) 
meV. 



gap function: A = 3.0(1. 43 cos 20 -0.43 cos 60) meV. The 
solid lines are the numerically calculated curves using Z 
= 3.0, T = 0.65 meV, a = for the spectrum along (100) 
direction, and a — it /A for the spectrum along (110) di- 
rection. One can see that, for the tunneling spectrum 
along (100) direction, the calculated curve coincides with 
the data at high bias voltages but significant deviations 
occur at low bias voltages. For tunneling spectrum along 
(110) direction, ZBCP is clearly seen in the calculated 
curve, in sharp contrast to the data. Therefore, the tun- 
neling spectra cannot be explained by d-wave gap sym- 
metry. 

Finally, the most powerful way to distinguish between 
any d-wave and anisotropic s-wave gap symmetries is to 
study the response of a superconductor to nonmagnetic 
impurities or disorder. The nonmagnetic impurity pair- 
breaking effect is both bulk- and phase-sensitive. This is 
because the rate of T c suppression by nonmagnetic im- 



FIG. 3: Comparison of the tunneling spectra of NCCO with 
the calculated curves based on a nonmonotonic d-wave gap 
function: A = 3.0(1.43 cos 29 - 0.43 cos 60) meV. The solid 
lines are the numerically calculated curve using Z = 3.0, V = 
0.65 meV, and a — for the spectrum along (100) direction 
and a = 7r/4 for the spectrum along (110) direction, respec- 
tively. 



purities or defects in a two-dimensional superconductor 
|29| is determined by the value of the Fermi surface (FS) 
average < A(k) >fs, which depends sensitively on the 
phase of the gap function. More specifically, the rate 
is proportional to a parameter x = 1 — (< A(k) >fs 
) 2 /< A 2 (fc) >fs- It is clear that \ — f° r isotropic s- 
wave superconductors while x = 1 f° r d-wave and g-wave 
superconductors. For the anisotropic s-wave gap: A = 
3.59(1 - 0.11 cos40) meV, x = 0-006. An equation to de- 
scribe the pair-breaking effect by nonmagnetic impurities 
(or defects) is given by [29[ 



0.122(M};)V 

T~r. 



)-*(^)], (3) 



where fifl* is the renormalized plasma energy [29j, |30( 
in units of eV, p r is the residual resistivity in units of 
/if2cm, and "F is the digamma function. The renormalized 
plasma energy can be independently determined from op- 
tical conductivity. Optical data of Pri.ssCeo.isCuO^j, 
indicate fiSl* = 1.64 eV (Ref. 311). We will use this value 
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of hfl* to calculate T c as a function of the residual resis- 
tivity in terms of both d-wave and an anisotropic s-wave 
gap function inferred from the Raman spectra above. 
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FIG. 4: T c as a function of residual resistivity in optimally 
doped Ndi.85Ceo.i5Cu04-i/. The data are from Refs. [Til. Il7i. 
l32l ]. The solid line is numerically calculated curve in terms of 
any d-wave gap and the dotted line is the calculated curve for 
an s-wave gap function proportional to (1 — 0.11 cos 49). 

Figure 4 shows T c as a function of residual resistivity 
in optimally doped Ndi.sBCen.iRCuCu-,,. The data are 
from Refs 
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32J. The solid line is the numerically 
calculated curve in terms of any d-wave gap and the dot- 
ted line is the calculated curve for an s-wave gap function 
proportional to (1 — 0.11COS46 1 ). For the d-wave gap sym- 
metry, the parameter-free calculation (solid line) shows 
that T c will be suppressed to zero at a very small residual 
resistivity of 10.7 /zficm while the measured T c is nearly 
independent of the residual resistivity. For the s-wave 
gap proportional to (1 — 0.11 cos 49) or to (1 — 0. 17cos46>), 
the calculated T c is nearly independent of the residual re- 
sistivity, in agreement with the data. Therefore, the data 
in Fig. 4 rule out any d-wave gap symmetry and unam- 
biguously point towards nodeless s-wave gap symmetry. 

In summary, the bulk-sensitive data of Raman scat- 
tering, optical conductivity, magnetic penetration depth, 
directional point-contact tunneling spectra, and nonmag- 
netic pair-breaking effect in optimally electron-doped 
Ndi.ssCeo.isCuC^-j, unambiguously support a nearly 
isotropic s-wave gap. The s-wave gap symmetry is consis- 
tent with the earlier [33| and recent [ij] conclusion that 
high-temperature superconductivity in electron-doped 
cuprates is mainly caused by electron-phonon coupling. 

* Correspondence should be addressed to 
gzhao2@calstatela.edu 
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